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It has been recently shown that meta-stable domain walls exist in high-density QCD (/i ^ 0) as 
well as in QCD with large number of colors (N c — > oo), with the lifetime being exponentially long 
in both cases. Such metastable domain walls may exist in our world as well, especially in hot 
hadronic matter with temperature close to critical. In this paper we discuss what happens if a 
bubble made of such wall is created in heavy ion collisions, in the mixed phase between QGP and 
hadronic matter. We show it will further be expanded to larger volume ~ 20/m 3 by the pion 
pressure, before it disappears, either by puncture or contraction. Both scenarios leave distinctive 
experimental signatures of such events, negatively affecting the interference correlations between the 
outgoing pions. 



I. INTRODUCTION 

Domain walls are common in field theory, the simplest 
in the family of topological objects. If they are configu- 
rations of fields interpolating between two distinct vacua, 
the domain wall obviously cannot decay. However even 
in theories with a single vacuum domain-wall configu- 
rations may nevertheless exist. Such domain walls are 
solitons which usually require a discrete-type of symme- 
try. In QCD such a symmetry indeed exists and corre- 
sponds to the discrete rotations of the so-called 6 angle 
9 9 + 2nn. 

The 9 parameter is not a dynamical field in QCD J 
and therefore it can not make a domain wall by itself. 
However, the 9 angle appears in the low energy descrip- 
tion of QCD together with the physical rj'(x) field in 
a very special combination (9 — i log DetU) where the 
unitary matrix U describes the pseudo-goldstone fields 
and the singlet rf{x) field. Therefore, the domain walls 
which may exist due to the discrete symmetry discussed 
above, can be realized by the physical rf{x) field. If 
the rj' field is the light field, all calculations are under 
theoretical control and one can argue that the life time 
of the domain wall is parametrically large in large N c 
limit Tuf e ~ exp(A 2 ), pi], as well as in large /j, limit 
Tuf e ~ exp(^ 2+b ),6 = yiV c - |iV/ {§ if one considers 
the high density QCD with large \i. 

In our world with iV c = 3 the rf meson is not a light 
particle, although it may become lighter for sufficiently 
hot hadronic matter as a result of partial restoration of 
U(1)a symmetry around the QCD phase transition, as 
argued in (3|J]. At the moment, we do not know what 
exactly happens with this excitation at such conditions, 
there is no theoretical control or sufficient lattice data. 
However, as was argued in Rcf. H , the domain wall in the 



theory with N c = 3 can still be classically stable. In re- 
ality (at quantum level), it is of course only metastable, 
but with a life time much longer than a typical QCD 
time scale, ~ ^-qcd- The argument of Ref. [j] is based 
on the observation that one should not naively compare 
rrirfi mass with a typical hadronic scale which is the same 
order of lGeV. Instead, one should compare the vacuum 
energy density due to the gluon degrees of freedom (it 
can be explicitly expressed in terms of the gluon con- 
densate, E = (ff^G 2 ) ~ TV 2 ) with the corresponding 
contribution due to the rf excitation when the relevant 
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Exactly this inequality prevents the domain walls from 
the classically allowed fast decay as discussed in [[I). In 
what follows we assume that rf domain wall is classically 
stable object, and therefore, it decays through the quan- 
tum tunneling process with exponentially large lifetime^ 
which is longer than any other time scales existing in the 
heavy ion collisions. 

The main point of this letter is the observation that if 
such domain walls indeed exist in QCD, they can be pro- 
duced and studied in heavy ion collisions. At high colli- 
sion energies (SPS, RHIC) the excited matter is assumed 
to be produced in the QGP phase, and then cools down, 
spending significant time ( ~ 5/m/c) in the so called 
mixed phase. Small bubbles made of domain walls can be 
produced by thermal fluctuations at this stage: currently 
we are not able to provide any quantitative estimates of 
the probability of this to happen, as it would require an 
understanding of the "out of equilibrium" physics. How- 
ever, one may not expect a huge suppression for the pro- 



*We do not discuss axions in this paper. 



t Recall familiar soap film bubbles: those are metastable as 
well, but they do exist for minutes with sizes of many cm, 
which is very impressive if expressed in terms microscopic 
(atomic) units. 
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duction of a ~ fewfm size bubble due to the general ar- 
guments that all dimensional parameters of the problem 
have one and the same ~ 1/m scale during this period. 
In fact the domain wall tension a ifO] : 



2V2N C 



(200AiW) 3 ~ (1/m)" 



(1) 



is not large enough to upset energetics at T = T c rs 
160 MeV. Our main observation is that after being pro- 
duced, a small (1/m in size) bubble would grow sub- 
stantially in size, driven by the pion pressure, spending 
a significant extra time (~ 5)/m/c as a relatively large 
macroscopic object, before releasing all its content and 
disappear. 

How can these large bubbles be observed if produced? 
Our suggestion is to monitor (on event-by-event basis) 
the strength of the two-pion Bose-Einstein correlation 
function, A*(p), see below, which it turns out to be very 
sensitive to any macroscopically large object with long 
enough life time. All pions which are eventually emitted 
from such an object will be completely incoherent with 
the rest of pions. 

Closing the Introduction we should mention that re- 
lated but different macroscopically large configurations 
were also discussed in Ref. [g) in the context of the decay 
of the metastable vacua possibly created in heavy ion col- 
lisions. The difference with this work is quite significant: 
we do not consider a metastable vacuum, but rather a 
metastable wall making a bubble. 



II. BUBBLE DYNAMICS 



Effective Lagrangian for bubble motion is 



L 



4tt<t R 2 (t) 



Air 

R 2 {t) - 4iraR 2 (t) + —R 3 {t)P v (2) 



where a is the wall surface tension ([I]) and R(t) is bub- 
ble radii. The last term describes the effective pressure 
induced by pions, which competes with the surface ten- 
sion and try to expand the bubble. It appears because 
pions are scattered back, from the domain wall into the 
bubble. The probability of this to happen P(k r ), the re- 
flection coefficient of 7r meson with momentum k r off a 
domain wall, enters the pion pressure 



P, 



2k r (k)P(k r ) 



d 3 k 
(2tt) 3 



(3) 



(P(k r ) will be estimated in the next section.) We assume 
the bubble is a spherically symmetric object such that 
we keep only the r component, v r is the velocity of the 
pions; n^ik) is the pion density, to be discussed below. 
The equation of motion is then 



R{t) --^-R + T' 



(4) 



which has perfect physical meaning as the last term de- 
scribes the acceleration of the bubble's surface. Indeed, 
Ptt can be interpreted as the force/area applied to the 
bubble, and a can be interpreted as mass/area of the 
bubble. 

The bubble radius R(t) is one time-dependent variable: 
the second one is the temperature of the hadronic gas in- 
side the bubble T(t). (We assume that pions rescattered 
back from the bubble wall are quickly equilibrated with 
the rest of pions inside it.) We therefore need the second 
equation, which we get simply from the energy conserva- 
tion. One can identify the following terms in the energy 
change of the bubble 



where 



Eqas P^leakaqe Email 



■ d ^R 3 {t) . 
E gas - -{c SB T 



and csb = 7r 2 /10 for massless pion Bose gasr 



Eleakage = ^R 2 / k r V r n n (k)(l - P{k r )) 



d 3 k 

2(tt) 3 



wall 



8tt<jRR 



(5) 



(6) 



(7) 



(8) 



After the replacement (v r k r ) — > in eq.(0) ( brackets 
imply the averaging over ensemble with temperature T), 
the corresponding equation finally gets the form 



AT R 
3T + R 



P* 



2aR 



1 



3R 2c SB T 4 R R 2 c SB T 4 



(9) 



If the rhs is set to zero, the bubble evolution gives just 
T 4 R 3 = const, as energy conservation would require. 
The negative sign of the rhs tends to compress the bub- 
ble. 

Our main observation is that after being produced, 
small bubble would be forced to grow substantially, pro- 
vided some conditions (to be specified below) are met. 
This expansion may happen due to two different reasons. 

The first, if a bubble happen to contain a signifi- 
cant portion of its volume by QGP, it will expand sim- 
ply because while it is transformed into hadronic matter 
(approximated by a pion gas) it occupies much larger 
volume. The simplest (but naive) estimate of such ex- 
pansion ratio R^ xp is obtained if both the QGP and 
the pion gas are treated as ideal massless gases; then 
it is just the ratio of degrees of freedom (DOF) of QGP 



(cqgp = 47.5, N c 



N 



f 



3) to that of the pion gas 



(^ = 3). If so, the volume is expected to grow by an 



■"■It is 9/-7T 2 in Boltzmann approximation we will use for 
simplicity. 
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order of magnitude. In reality, we have the resonance 
gas, not the pion one, with more degrees of freedom, 
and partons are not really massless. So R exp is not 
cqgp/ctt ~ 16, but somewhat smaller. 

The second mechanism driving expansion works in 
hadronic phase. We will show in the next section that 
the thermal pions are rather effectively reflected by the 
domain wall, so that they are effectively trapped in- 
side the bubble. It leads to another expansion ratio 
B% cp = s(T c )/s(T m i n ) where the entropy densities at two 
temperatures are included. The latter one, T m , n , should 
correspond to final mechanically equilibrium bubble sat- 
isfying the following condition 

P {T m m) = a/R f (10) 

The larger the final bubble radius Rf, the lower this pres- 
sure. The lifetime of this equilibrium is determined ei- 
ther by partial leakage of pions, or by the bubble lifetime, 
whatever is shorter. 
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FIG. 1. An example of the solution of equations of motion 
of the bubble dynamics, described in the text. The solid line is 
the dependence of the bubble radius R(t), in fm, versus time 
t, in fm. We assume R(t = 0) = 3/m which corresponds to 
the initially small bubble about 1.2/m expanded due to the 
first mechanism, 1.2/m ■ ( CQGP )V 3 ~ 3/m up to the rela- 
tively large size R(t = 0) ~ 3/m. The dotted line represents 
the corresponding evolution of the internal temperature T (in 
fm- 1 ). 

Equations (|J) and (||) derived above, describe the dy- 
namics of bubbles in this hadronic phase and correspond 
to the second stage of expansion determined by R e ^ p '. 
We solved equations @) and (||) numerically to demon- 
strate the effect of this expansion, see Fig.l. Let us only 
comment on T evolution. It starts at the critical value 
T c = 165 MeV — .83 fm -1 and then cools down as the 



bubble expands. However, the equations indicate a sec- 
ondary heating as the bubble walls collapses. We will not 
dwell on it, but note that quite similar phenomena are 
known elsewhere in physics, e.g. in the so called sono- 
luminescence: so it may be not just an artifact of the 
equation solution. 

III. PION SCATTERING ON A DOMAIN WALL 

In the discussions presented above we introduced in 
our formula (||) one essential parameter, the reflection 
probability P{k z ) which is, by definition, the probability 
of reflection of 7r meson off a domain wall. The behavior 
of this parameter has not been specified yet, we shall 
estimate P{k z ) now. 

First of all, we start from the low energy effective La- 
grangian when 7r and rf fields are described by the uni- 
tary matrix in the simplified version of the theory when 
Nf = 2, m u = m d : 

^(^+^,^ = 1,, (11) 

where cr° are the Pauli matrices, 7r a is the triplet, and 
/„. ~ /„/ ~ HZMeV . In terms of U the low energy 
effective Lagrangian is given by || : 

L = f J-Tr{d^d^U) + \MTr{U + Eft) + 
o 2 

„ .HogDetU — 6. , „. 

£cos( — ), (12) 

where all dimensional parameters in this Lagrangian are 
expressed in terms of QCD vacuum condensates, and are 
well-known: M = mq\{'ij}ip)\; the constant E is related to 
the gluon condensate E — (^f^G 2 ). The first two terms 
describe the standard expression for the effective chiral 
Lagrangian; the last term describes the 7/ field. This 
term is not uniquely fixed by the symmetry; however the 
only important element for the following discussions is 
the manifest 2tt periodicity for the r/' term as well as 
an appearance of the scale E which remains finite in the 
chiral limit. A specific cos form for this term is not essen- 
tial; however it will be used for all numerical estimates 
presented below. 

As we discussed in the Introduction the theory jl^ ) 
supports the metastable domain walls. We refer the 
reader to the original paper Q for the discussions of the 
specific properties of the domain wall; the only what we 
need to know here about the QCD domain walls is their 
following general properties: 

a). The QCD domain wall is described by two dimen- 
sionless phases, (ps(z) describing the isotopical "singlet", 
and <)>t{z) isotopical "triplet" fields. These fields corre- 
spond to the dynamical rj' (singlet) and pion ir° (triplet) 
fields defined in (pd|); they depend only on one variable 
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z. Both dimensionless fields <fis(z), 4>t(z) interpolate be- 
tween zero and 2tt when z varies from — oo to +00 and can 
be expressed in terms of the original dimensional fields 
rj', 7T° as follows. 



5 S,T 



4>u ± <t>d, v' 



cf> S (z), 7T° = ^=MZ)~ (13) 



2V2 



2V2 



b) . Metastable domain walls of a minimal energy cor- 
respond to the transitions from (4>m4>d)\z=-oo — (0,0) 
to (cj> u , 4> d ) I z=+OQ = (2tt, 0). In the limit m u = m d the 
transition to (<fi u , 4>d) = (0, 27r) have the same energy and 
there is a degeneracy. In reality md > m u , and the tran- 
sition to ((j) u ,4>d) = (27r,0) is the only stable transition. 

c) . Domain wall has a sandwich- like structure: <fts sub- 
stantially varies on the scale z ~ m^, 1 while <pT varies 
on considerably larger scale z ~ m' 1 ^> m~} where <ps 
is already close to its vacuum values 0, 2tt. Therefore, 
one can say that the 77' transition is sandwiched in the 
pion transition. In spite of the fact that </>t is much 
wider than (f>s the main contribution to the wall tension 
(|l|) comes from the <ps transition. Anti-soliton corre- 
sponds to the transition from (<fi u , 4>d)\z=-oo = (0,0) to 
(4>u,MU=+oo = (— 2tt,0). 

Our next step is to consider small oscillations in the 
pion and rj' fields about the static QCD domain wall so- 
lution^): 

7r(aVt)-> (^r 1 {x l ,),Tr2(x ll ),^=(f> T (z) +tt 3 (x^ ; 



2V2 



where 4>t{z) and (j)s(z) are solutions known numerically 
and qualitatively described above. It is clear that a 
particle scattering by a planar domain wall reduces to a 
one-dimensional scattering problem. To further simplify 
things we neglect scattering of rj particles which can not 
play an important role due to its larger mass. 

First, let us consider the simple case of scattering of 
the neutral tto meson off the domain wall. With the plane 
wave ansatz, ttz(x^) — ^3(2) exp(— iuit + ik x x + ik y y), 
the field equation for 773(2) which follows from the chiral 
Lagrangian (113), can be written as 



(15) 



+k 2 z -U^ (z) ) ,r.C) 0. 
k 2 z = u?-k 2 x -k 2 y - ml, 



where the effective potential U w o (z) for the problem is 
expressed in terms of the domain wall profile functions 
4>t(z) and <ps(z): 



U 7! o{z) = — m w 1 — cos 



4>t(z) 4>s{z) 



cos • 



(16) 



Equations (|15|), ( |16| ) can not be solved analytically, how- 
ever, a qualitative behavior of the reflection coefficient 



P{k z ) can be easily understood from the following ar- 
guments. If the particle wavelength is much smaller 
than the thickness of the wall, k z ^> m w , the reflection 
coefficient is exponentially small according to the stan- 
dard semiclassical arguments. In the opposite, the long- 
wavelength limit k z — ► 0, the wall potential ( |l6|) can be 
adequately approximated by a 8(z) function, 



U^o(z) ~ —aS(z), 



a 



U„o(z)dz ~ m w . (17) 



The scattering problem with a 8{z) function potential is 
easily solved; the reflection coefficient is 



P{k z ) = 



{2k z f 



(18) 



As expected, the reflection coefficient does not depend 
on the sign of the potential U v o(z), and P(k z ) — ► 1 for 
small k z — > 0. 

Our next task is an analysis of a similar problem for the 
charged 7r + ,7r~ components of 7? field. Some technical 
complications arise here due to the fact that the kinetic 
term ~ Tr(d fl W d^U) in the low energy Lagrangian ( |12) 
is not reduced to the canonical kinetic terms ~ l/2(d u n) 2 
for the individual components when the expansion in (FL3) 
is made about a nontrivial classical configuration. In- 
stead, the kinetic term is given by 



-Trid^d^U) = (d,6) 2 + (8^ 



+ sm 2 0{d ll n) 2 , (19) 



where, in order to simplify formula (|l9|), we introduce 
new variables 9{x), 4>{x), n(x) expressed in terms of the 
original fields 7?, 77' dTll) as follows: 



= V 7? 2 ; n = — ; 



_ V2\tt\ 



V2rf 



(20) 



After some algebra we arrive to the following expres- 
sion (analogous to ( |l5| , |l6| )) describing the scattering of 
charged tt± , 1:2 components off the wall. 



— + k 2 z -U n ±(z))7T t (z) = 0, 1 = 1,2 



k\ = uj 2 - k 2 x - k 2 y - ml, 



(21) 



where the effective potential U„±{z) for the charged 771,2 
components is expressed in terms of the same profile func- 
tions </>t(z) and (f>s( z ) as follows 

U^{z) =ml H Z)COS ^ -l) +SUM*l (22) 



where 5U 7T ± (z) is a quite complicated operator defined 
in the Appendix and numerically will be ignored for the 
qualitative discussions which follow. Note, that in all for- 
mulae presented above the profile functions (f>s{z), 4>t(z) 
describing the wall are defined in the region — 00 < z < 
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when < (ps(z),(j>r{z) < it. For the positive z, by sym- 
metry, one should replace 4>s.t —> 2-7T — (f>s.T as explained 
in p]]. Qualitative properties of the potential U v ±(z) are 
the same as U n o(z) discussed earlier (16), namely, a long- 
wavelength particle is reflected from the wall with a very 
high probability, while in the short-wavelength limit the 
reflection probability is nearly zero. The reflection co- 
efficient is slightly different for the charged and neutral 
pions; however, in what follows we neglect that difference 
and approximate the reflection coefficient as follows 



P(k z ) ~ e"tf 



(23) 



This is our final expression which was used in the previ- 
ous section, see eqs.(^), (|7|), for the numerical estimates 
presented on Figl. 



IV. EXPERIMENTAL OBSERVABLE 

Our main idea of the experimental observation of bub- 
bles is to make use of the intensity interferometry of pi- 
ons, due to their Bose-Einstein statistics. This method, 
also known as HBT (Hanbury-Brown-Twiss) interferom- 
etry, has been originally introduced for measuring the 
angular diameters of stars. For pions it has been first 
observed in pp annihilation, and explained in the famous 
paper by Goldhaber et al j7j. Since early 70's it has been 
used to extract source sizes of hadronic fireballs, see e.g. 
one of the early papers on the subject by one of us 
Recently it has been argued M that one can use pion 
interferometry as a sensitive tool to detect possible in- 
crease of the ?y' production in heavy ion collisions. We 
now extended the same reasoning for the observation of 
the bubble production. 

If a bubble is produced, it exists for some lifetime and 
then decay. It can either happen due to (i) puncture of 
the wall, or (ii) simple contraction, as discussed above. 
Both scenarios lead to similar observable signature. 

In the case of puncture, the bubble walls decays into its 
underlying fields, the 77'. Since it happens after a freeze- 
out time for most pions not belonging to the bubble, most 
of these 77' would not be re-absorbed, and decay normally, 
most often into 5 pions. One may argue (direct following 
ref. H) that there exist an important observable distinc- 
tion between those pions and the rest of them produced 
from the fireball. Large lifetime of rf ,77 make their prod- 
ucts incoherent to other pions, reducing the HBT peak 
in two-pion spectra. It happens, because the inverse of 
the 77' lifetime is much smaller than experimental energy 
resolution of the detectors. 

The second (and, as our estimates suggest, see Sec- 
tion II, more probable ) mechanism for the bubble de- 
cay is described above: it is due to a eventual collapse 
of the bubble due to pion leakage. In this case bubble 
surface contracts and large population of n' is not ex- 
pected. However this processes takes long enough time, 



~ 5 fm/c, so the bubble itself plays the role of long- 
lived object. Again, the inverse of this time is likely to 
be below the detector resolution. As a result, most of 
the pions trapped in these quasi-stationary bubbles and 
released later become incoherent with others pions. For 
this reason the pions from bubbles lead to the same effect 
of not producing a HBT peak in two-pion spectra. 

The strength of the HBT correlation is characterized 
by the effective intercept parameter A* (j>) defined as fol- 
lows § 



A.(p) 



N, 



direct 



(p) 



direct (p) + Ndelayedip) 



(24) 



where N direct (N de i ayed (p)) is the one- particle- invariant 
momentum distribution of the "core" direct ( "halo" or 
delayed ) pions. As was discussed in 0, a substantial 
increase in the rf production will result in appearance 
of a hole in the low pt region of the effective intercept 
parameter A*(p) centered around pt — !38MeV which 
represents the average pt of the pions coming from 77' 
decay. If a bubble is punctured, we also expect the same 
signature, due to an increase number of 77' to be produced 
with low p t . 

However, if the second mechanism of decay (slow infla- 
tion and deflation of the bubble) is prevailed, the effect of 
decreasing A*(p) is different. This is because the emission 
of pions from the bubble happen through the walls, and 
low energy pions are expected to be nearly completely 
trapped (note, the reflection coefficient is close to one for 
the cold pions and zero for the hot ones, see previous 
Section III). So we expect to see a decreasing A»(p) at 
larger p instead. 

To make a numerical estimate of the effect, we note 
that the parameter A* equals to 1 for completely coher- 
ent pions and reduced to about 0.5 in usual experimental 
conditions^} If bubbles are produced, the intercept A» 
would be additionally reduced by the factor (1 — f bubble) 2 
as follows from ( pjh . Here f bubble is the fraction of pions 
coming from the bubbles, this parameter can be eas- 
ily estimated as follows. The bubble energy is order of 
Ebubbie — 4:nR 2 o- ~ 60GeV, where we use R ~ 5/m and 
a ~ lfmT 3 (Q). If all the energy accumulated in the wall 
of the bubble will go to the production of the 77' mesons 
(which will result in additional ~ 3O77' mesons per event 
) one should expect a 100 or so of "incoherent" pions 
to be produced from the bubble. In the second scenario 
(inflation/deflation) for bubble decay the effect would be 
proportional to i? 3 rather than R 2 and because all the pi- 
ons from the bubble incoherent, therefore their number 
could be similar or even larger. 

Naively this number represents relatively small fraction 
of the total number of pions in each given event. How- 



5 There are delayed pions due to "natural cocktail" of long 
lived resonances. 
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ever, the bubble is local in rapidity. Even at mid-rapidity 
at RHIC, with total pion multiplicity about 1000 per unit 
rapidity, the products of "naturally occurred" long lived 
resonances make about 300 of those. Additional 100 in- 
coherent pions expected from the bubble makes it 400, 
with the average A changing from .7 2 = 0.5 to .6 2 = 0.36, 
not a non- negligible effect. 

Therefore, we propose to look at the event-by-event 
fluctuations of parameter A*, hunting for the tail of the 
distribution toward its values smaller then the average, 
preferably at low p t . An unusually long tail may indi- 
cate the bubble formation. The thermodynamical fluc- 
tuations in the particle composition are expected to be 
very small 0(l/y/N), effect, obviously unable to produce 
any long tail by itself. 
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APPENDIX A. 



Here is the definition of SU^± (z) which appears in 
eq.(B2) in terms of the profile functions 4>s(z), (j>x(z): 



2 sin 



dz 



SU T ±(z)n= (25) 



/ 



dz dz \4>t dz 



with f(z) defined as 
f = l- 



2 . Mz) 
— sin 

vt 2 



(26) 
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